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04 ' We construct lattice parafermions - local products of order and disorder operators - in nearest- 

neighbour Zjv models on regular isotropic planar lattices, and show that they are discretely holo- 
^ . morphic, that is they satisfy discrete Cauchy-Riemann equations, precisely at the critical Fateev- 

Zamolodchikov (FZ) integrable points. We generalise our analysis to models with anisotropic inter- 
actions, showing that, as long as the lattice is correctly embedded in the plane, such discretely holo- 
, morphic parafermions exist for particular values of the couplings which we identify as the anisotropic 

' FZ points. These results extend to more general inhomogeneous lattice models as long as the cov- 

ering lattice admits a rhombic embedding in the plane. 



I. INTRODUCTION 

, Zat models, as the simplest statistical mechanics models with discrete global symmetries, have been studied for many 
years. They are interesting, both as lattice models and in the continuum limit as quantum field theories, because 
they exhibit semi-locality properties, such as fractional spin ^Jj, which arise in various other domains of physics. 
The simplest examples of these models are the well-known Ising model (A^ — 2), three-state Potts model {N — 3), 
] and Ashkin- Teller model (A^ = 4), all exactly solvable to some extent. The phase diagram of the nearest neighbour 
^ . Zjv model in two dimensions for A" = 5, 6, 7 was elucidated in 0, [^, but in the general case it is still not known. 

Although these lattice models for A^ > 4 have complicated phase diagrams with several critical surfaces, Fateev and 
^ — , Zamolodchikov [3] showed that there are some points (hereinafter referred to as FZ points) in the critical surface 

■ at which these models are solvable, in the sense that they satisfy generalised star-triangle relations. Moreover there 
is strong evidence that the scaling limit at these points corresponds to Zamolodchikov and Fateev's parafermionic 
conformal field theory f^. In this article we give further evidence for this connection. In particular we find lattice 

1^ I candidates for the holomorphic parafermions of the continuum model, and show that these satisfy a discrete version 

■ of the Cauchy-Riemann equations. We also extend this result to models with anisotropic interactions. We show that 
as long as the lattice is correctly embedded in the plane, we can again construct discretely holomorphic parafermions. 
We may use this to conjecture the location of the FZ critical surface of these anisotropic Z^v lattice models. In fact, 
our methods can be extended to any non- uniform lattice of the Baxter type thanks to a result of Kenyon and 
Schlenkcr 7] about embedding such lattices in the plane. 

Identifying discretely holomorphic objects in lattice models is an important step in showing that the scaling limit of 
suitably defined curves in these models is described by Schramm-Loewner evolution (SLE) [1, Q. Suitable candidates 
for SLE curves in Zjv models have been suggested by Santachiara [l^l (see also Gamsa and Cardy [llj for the case 
N = ?>). We expect that our results will be the first step in showing that this is indeed the case. 

This paper is organised as follows. In the second section we briefly review some related properties of Z^r models, 
and the corresponding parafermionic conformal field theories. We define the disorder variables on the isotropic square 
lattice, similarly to the case of the Z n clock models li] and the Potts model . We then define lattice parafermions 
as products of neighbouring order and disorder variables, with a suitable phase factor. Using the Boltzmann weights 
of the model at the FZ point, we then show that there are local linear relations between parafermions, which are 
equivalent to discrete holomorphicity. This is our first main result. 

In the third section we consider lattices with anisotropic interactions and show that there are special points, which 
we identify with anisotropic FZ points, at which there is a certain embedding of the lattice on the plane for which 
discrete holomorphicity is again recovered. We further extend this to the case of Baxter lattices @ , where we show 
that the existence of a rhombic embedding 0] for such lattices allows us simply to generalise our results. It should 
be mentioned that the importance of rhombic embeddings (and the related isoradial embeddings) was central to the 
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work of DufRn [T^ , Mercat [T3| and Bobenko et al. [Tsl in the general theory of discrete holomorphy. Similar ideas 
have been used in other lattice models by Kenyon [161] and Bazhanov et al. [17| . 

II. HOLOMORPHIC PARAFERMIONS IN Z(N) MODEL 

The general Z^r model on an arbitrary graph can be defined by associating to every node of the graph a variable 
Sr which takes values in the set lj'', g = 0, 1, TV — 1, where uj = e^'^*/^. Each configuration occurs with probability 

Z-'\{W{ry), (1) 

(rr') 

where the product is over all the edges (rr') of the graph, and Z is the partition function. The weights W{r, r') take 
the general form 



W{ry)=Y,xt'\s^s:,)\ (2) 

fc=0 



where * denotes complex conjugation. Assuming x''^^ 7^ 0, we can always set it equal to unity. Reality of the weights 
implies that x^^^ ^ = a;^^]., which means that one can describe them by just [N /2] real parameters. The weights are 
invariant under the global group Z^r, Sr Lo^Sr, as well as charge conjugation C : Sr ^ s*. 

In what follows, except for the last section, we consider the graph to be a regular lattice, with translationally 

invariant couplings x^^^ \ which may however depend on the lattice orientation of the edge (rr'). 

For N = 2 and 3 it is easy to see that the model is equivalent to the Ising model and three state Potts model 
respectively, whose critical points correspond to well understood minimal conformal field theories with c < 1. The 
higher values of N are related to more complicated theories with rich phase diagrams [Tsj. Nearest-neighbour Zn 
models exhibit Kramers- Wannier duality symmetry, which is important in understanding their phase diagrams. If 
all of the Xk with A; > 1 have the same value, we have the well-known iV-state Potts model, which has a single 
first-order transition for N > 4. In general for TV > 3 the theory has a critical surface on which the exponents may 
vary continuously. It has been conjectured that there are n < points in the phase diagram which correspond to 
special kind of conformal field theories called Fateev-Zamolodchikov (FZ) parafermionic models of the first kind 
This has been verified by numerical simulations ^ISj . Our lattice holomorphic operators, described below, are another 
evidence for this conjecture. These conformal field theories have the central charge 

2(TV-1) 
N + 2 ' 

the cases iV = 2, 3 coinciding with the central charge oi p = 2, 5 minimal CFTs. Kramers- Wannier symmetry 
manifests itself in — 1 order parameter fields am and N — 1 dual disorder fields /im, which are conjectured to be the 
continuum version of the order and disorder variables on the original lattice. Zamolodchikov and Fateev, using only 
the Zjv symmetry, showed that there also exist holomorphic operators (that is, operators whose correlation functions 
are analytic functions except at coincident points) with the following dimensions (equal to their conformal spin) 

m(N — m) , ,r , /.X 

TV ' m = l,...,Ar-l. (4) 

We now give our construction of the lattice analogs of these objects. We first consider the homogenous isotropic 
Zjv model on the square lattice. 

In the first step let define the disorder operators on the sites of dual lattice, which, in this case, is again another 
square lattice whose vertices are at the centers f of the faces of the original lattice. The insertion of a disorder operator 
fifm corresponds to modifying the weights so that the order operator Sr has monodromy Sr — > uj~"^Sr on taking the 
point r in a closed circuit around f. This is equivalent to introducing a path, or string, on the sites of the dual lattice 
from r to infinity (or some other point on the boundary), such that the weights on edges (rr') intersected by the 
string are modified by the substitution SrS*, — > SrW~™s*/. 

Thus the disorder operator has the following form for the general Zjv model 

{rr') intersected by string ^k=o XkySrS^., ) 
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FIG. 1: An elementary square of the covering lattice (solid lines - rotated by 45° to correspond to the conventions in the text.) 
The opposite pairs of vertices are associated with order (s) and disorder (fi) operators respectively. The parafermions ^ are 
associated with the edges. Discrete holomorphicity at the FZ points means that the contour sum of the parafermions around 
each face vanishes. Also shown is part of the original square lattice (dashed lines) and the string (dotted line) attached to the 
disorder operators. 



It is not difficult to sec that the disorder variables have the same C-symmetry similar to the spin variables, 
fJ-rn = fJ-N-rm and that, up to a gauge transformation, the definition ^ is path- independent. 

Consider the square which the vertices are made by the two neighbouring spin variables si, S2 and the two neigh- 
bouring disorder variables /i];,/i2 (Figl.H]). By taking the string from 1 to run along the dual edge (12), from ^ we 
can write the following relation 

Mim = ^N~l : ^^27n ■ (6) 

We will need the exact value of the Xk at the FZ point, which for the isotropic square lattice have the following 
compact form Q 



fc-i 



•^ck 



j=0 ''"H AT 4n) 



For clarity, let us first focus on the simplest case N — 2. By multiplying both side of (O by the denominator and 
then by multiplying the result by si we find the following equation 

Sifll + XcS2til ^ Sifl2 ~ XcS2fJ-2 ■ (8) 

A similar equation can be found by exchanging si <-^- S2- By combining these two equations (multiplying the second 
equation by —i and adding this to the first) one finds 

— sifj-i — ie~~ S2fJ'i + S2M2 + ie~ sifj,2 = . (9) 

Note that for this to happen it is crucial that we use the critical value Xd — tan(7r/8) = \/2 ~ 1- (For the special case 
of the Ising model, when x ^ Xc the right hand side of ^ is proportional to the antiholomorphic fermion Q.) 

^ has the form of a discrete contour integral around each elementary square of the covering lattice (the union of 
the dual lattice with the original lattice) 

y IpeSZe — with Vrf — C ^ Srfif , (10) 

e 

where 9rf is the angle that the directed segment rr makes with the s-axes, with the convention — tt < 9rf < tt. In this 
case these angles are just —tt,—^,0, ^ but in the next section we will see more general cases. By changing i ^ —i 
one can find the similar equation for the discrete antiholomorphic fermions. These two linear equations are a lattice 
discretised form of the Cauchy-Riemann equations. 
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The calculation for = 3 is very similar. Multiplying both sides of equation ^ by the denominator, and then by 
si, S1S2 and S2 we can find three linear equations in the six variables si^i 2, SiS2Mi 2' *2/^i 2- -^Y eliminating S2/ij 2 
one again finds a discretely holomorphic equation of the form pop with now 

Iprf = e 3 Sr/if . (11) 

Once again, this works only at the FZ point given by ([7|). 

The case iV = 4 is very similar and we will discuss the results later. = 5 is more interesting. This model has two 
FZ points {xi,X2) = {xci,Xc2) and {xi.2,Xci) corresponding to the transformation which exchanges Xi and 

X2. For this case the equation ^ has the following form 

/i];„(l + XclSiS*2 + Xc2{siS*2Y + Xc2{siS*2)^ + Xcl(siS2)'') 

= A^2m(l + ^ ~XclSlS2 + e =~a;c2(siS2)^ + e~^~^ Xc2{siS*2)^ + e~^^ Xcl{siS*2Y) . (12) 

For the m — 1 case by multiplying the above equation by si, S2, S1S2, s\s'^ , Sis'2' finds the following equation 

/ 1 2*1 3 *2 I 4 *3 I *4\ 

PlI^{Si+ XclS-^S2+ Xc2SiS2 +Xc2SiS2 + XclS2 ) 

= /^2i(*i + e^^a^ciSiS2 + e"^a;c2S?S2^ + Xc2s\s*2 + e^^a;ciS2) , (13) 

and four other equations by permuting the terms si, S'{s2, Sis'^ , s\s'^ , S2- The next step is to eliminate the terms 
involving 3^82, Sis'^ , s\s^ from the above five equations. The calculation is rather cumbersome, so we resorted 
to using Mathematica. The result is as follows: for to = 1 it is possible to eliminate the above terms only at 
the first FZ point {xci,Xc2)^ in which case we obtain the holomorphic equation (jlOp with the holomorphic variable 

For TO = 2 we first multiply the equation p2|) by s\, s\s2t s\s'^ , s"^ , siSj"* to find 

/ 2 I 3 * , 4 *2 I *3 I *4\ 

lJil2{S-^+ XclS-^S2+ Xc2S-^S2 + Xc2S2 + XclSlS2 ) 

= At22('Si + e~^Xcls\s*2 + e'~^Xc2s{s*2 + e^^Xc2S*2 + XclSiS*2^) , (14) 

and four other equations by permuting the different terms. Again we should eliminate the terms in involving variables 
s\s2: s\s^ ^ 8182'^ . In this case we can find a holomorphic equation only at the second FZ point {xi,X2) — {xc2,Xci) 

with the holomorphic variable iprf = c 5 s^/Xf2- 

The calculation for to = 3 and m = 4 is similar and gives the antiholoniorphic parafermions. 

The generalization to higher values of N is straightforward. By simplifying equation ([6|) at the FZ points we 
conjecture that the following variables will be discretely holomorphic: 

^;:i = e-*^''"^-4>™, (15) 

for 1 < TO < [N/2]. The same expression with the argument of the exponential having the opposite sign is discretely 
antiholomorphic for [N/2] < m < N — 1. (This is just the complex conjugate of V'^~™-) We checked the conjecture 
for = 6. In this case we find all of the holomorphic operators at the one FZ point {xci,Xc2,Xc3)- These discretely 
holomorphic quantities are the candidates for the holomorphic parafermions of FZ parafermionic CFT in the continuum 
limit. Wc should emphasise that these operators arc holomorphic only at the FZ points and not elsewhere on the 
critical self -dual surfaces. For example, for = 4, ip^f is not holomorphic at the critical point of the four-state Potts 
model. However ip^f, which has conformal spin 1, is holomorphic all the way along the critical line 2xci + Xc2 = 1- 

We should also note another curious feature, which is most simply illustrated for the Ising case. Although we 
have shown that the linear relation ([9]) implies the discrete holomorphicity relation (fTO|) for the parafermion 4'rf = 

e 2 Sr^J^f, if we define another quantity iprf = e 2 s,./if:, then ^ also implies that 2_/e ^eO-^e = Oj that is ^prf 
is discretely antiholomorphic. This is quite general, and holds for other values of A^ and the more general lattices 
discussed in the next section: if tprf is discretely holomorphic, then ■0,.^ = e^^^^^'^^il^rf is discretely antiholomorphic. 

III. Zat model on other lattices 

In the previous section we found additional evidence that the Zamolodchikov-Fateev parafermionic CFT is a good 
candidate for the continuum limit of the FZ point of Zn models on the square lattice, by identifying discretely 
holomorphic parafermions on the lattice. Note this construction was possible only at the FZ points. One may therefore 
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try to invert the argument and locate the FZ point by requiring discrete holomorphicity. Since, by universahty, the 
same conformal field theory should describe the continuum limit of suitable Z^v models on other lattices, one would 
expect to be able to identify discretely holomorphic objects in this case, and thereby deduce the location of the FZ 
points. In fact we will show that this is possible in many cases, and that the critical weights of the models are just 
related to the geometry of the covering lattice when it is suitably embedded in the plane. Note that the notion of 
holomorphicity is not invariant under a general diffeomorphism of the plane, only under conformal transformations. 
Therefore we expect that the identification of discretely holomorphic objects will depend on choosing a particular 
embedding of the lattice into the plane, modulo conformal mappings. In what follows we study this problem for the 
anisotropic square, honeycomb and triangular lattices, and then the more general case of a 'Baxter lattice'. 

Let first investigate the square lattice with unequal weights a;|, x\ in the x and y directions respectively. In this 
case it is clear that in order to maintain invariance under lattice translations and reflection symmetry about the x 
and y axes the only transformations allowed are relative scalings of the x and y coordinates pL9]. In this case each 
elementary square of the covering lattice is deformed into a rhombus (Fig. [2]). 



1*/ 










IS, 



FIG. 2: Embedding of the anisotropic square lattice (dashed lines) in the plane. The faces of the covering lattice form rhombi 
with opening angles a (shown) and tt — a. 



Defining the parafermions by (jlSp . we can try to demand that the discrete contour integral around each elementary 
rhombus vanish as before. As an example we show the case of the Ising model in detail. Taking an arbitrary linear 
combination of ([5]) and its counterpart with si <-^- S2 we have 



y H y TS2^^l + ySl^.2 + S2A*2 = 16) 

ax\ — 1 axi — 1 1 — axi 

On the other hand, for the rhombus in Fig. [2] the contour sum is 

- Vii + e-*(^-")V2i + ^-22 + e^^V-iS = (17) 

where il^il = e^'^^^si/ij, -021 = e~'"/^S2Mi: "022 — ^'^^^2 ^^^^ V'lS — e*^'^~"^/^si^2- Comparing (fTB]) and pT|) we find 
that they are consistent only at a particular value of a and for x\ — tan(Qf/4). Similarly we find the contour sums 
around the rhombi aligned in the other direction can vanish only if = tan((7r — a)/ A). 

Note that the comparison of ()16|17p implies 3 complex equations. If instead of choosing the elementary face of the 
covering lattice to be a rhombus we take an arbitrary quadrilateral with edges 5ze, where ~ this would 

introduce 2 further complex unknowns (ratios of the 5ze) into (jl7p . besides the unknown a. Given x\, the 3 complex 
equations therefore determine a and the shape of the quadrilateral. We have already shown that a rhombus with a 
suitable opening angle satisfies all the equations. Therefore one would expect that they cannot be satisfied for any 
other quadrilateral with unequal edge lengths. This can be checked explicitly. 

For higher values of N we find that the contour sum around the rhombi vanish as long as the weights satisfy 

x^ - X Ja) = TT ^^"^^ + 

xl{a) = Xcfc(7r-a). (18) 
The result agrees with that found by Fateev and Zamolodchikov Q by imposing the the star-triangle relations. 
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Next we consider the isotropic honeycomb and triangular lattices, which are mutually dual with the vertices of the 
honeycomb lattice at the centers of the faces of the triangular lattice and vice versa. The elementary cells of the 
covering lattice give a regular rhombus tiling of the plane, with a — ^ for the honeycomb lattice and ct — for the 
triangular lattice. The mathematics then proceeds just as for the anisotropic square lattice above, with the result 
that one can find discretely holomorphic parafermions as long as the weights satisfy (jl8p with the appropriate values 
of a. This we then conjecture to correspond to the FZ point(s) for these lattices. 




FIG. 3: Embedding of one face of an regular anisotropic triangular lattice (dashed lines) in the plane. The faces of the covering 
lattice are rhombi with opening angles qi, 02 and Q3 = tt — ai — 02. 

A more interesting case is that of the homogeneous triangular lattice with unequal weights xl^\ x^^^ and xf^ in 
the three lattice directions. In this case, the question of where to locate the dual vertices is crucial. On the basis of 
the above observations, we adopt the following construction. Consider a (proper) embedding of the regular triangular 
lattice into the plane by some general linear transformation of the coordinates. This gives a regular tiling of the plane 
by triangles. For each triangular face with Z^v spins at the vertices, locate the dual vertex at the circumcenter , that 
is the point at which the three perpendicular bisectors of the edges meet, equidistant from the three vertices. This 
construction guarantees that adjacent pairs of vertices and dual vertices always form a rhombus. Each triangle is 
associated with three different rhombi with angles ai, a2, and = tt — ai — a2, see Fig. [31 If now we choose the 
weights so that 

xik = Xck{ai) , X2k = Xck{a2), X'ik = Xck{T^ - ai - a2) , (19) 

it follows from our general analysis that we then can identify discretely holomorphic parafermions. Transforming back 
to the orginal regular triangular lattice with unequal weights, we therefore conjecture that (|19p gives the general critical 
surface of the FZ point(s). Note that it is possible for one of the aj to be negative, if the largest angle of the deformed 
triangle is obtuse, corresponding to the interactions in one of the lattice directions being (weakly) antiferromagnetic. 
(However these points are still on the ferromagnetic critical surface, and have the Fateev-Zamolodchikov CFT Q as 
their conjectured scaling limit.) On transforming back to the regular lattice, the images of the dual vertices do not 
in general lie at the centers of the triangular faces (in fact they may lie outside the face). 

It is easy to check that these results are consistent with previously known results for the N — 2 Ising model and 
TV = 3 three states Potts model 

The FZ critical points for the honeycomb lattice with anisotropic couplings follow directly by duality from those of 
the triangular lattice, which corresponds to letting aj ir — aj in (|19p . 

We now indicate how our results extend to a more general kind of inhomogeneous lattice, called a Z-invariant or 
Baxter lattice This is a planar graph £ which is a union of AI simple (non-self-intersecting) curves crossing the 
complex plane from Xj — ioo to Xj -|-zc», where the values {xj} and {xj} (1 < J < M) are distinct, and with the further 
property that a given curve can intersect any of the others at most once (Fig. |4]). The faces of C are 2-colorable, 
and in general we can assume that its vertices are all of degree four (if not we deform the curves slightly so this is 
true.) Consider now the planar graph Q whose vertices are associated to each black face of £, and whose edges E pass 
through the vertices of C. We can define a Z^v model on Q with general weights x\ . The vertices of the dual lattice 
0* then correspond to the white faces of £. The vertices of the covering graph C (the union of the vertices of G and 
G*) each correspond to a face of £, irrespective of color. Note that all of the faces of C have degree 4. We are free to 
embed this lattice in the plane in any way we choose, just as for the anisotropic triangular lattice earlier. However, 
a remarkable theorem due to Kenyon and Schlenker Q states that (in the case considered here, where the original 
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FIG. 4: Part of a Baxter lattice. The faces of the graph L formed by the curved hnes are 2-colorable (not shown). Order 
variables s and disorder variables p are associated with alternately colored faces respectively. The covering lattice is shown as 
solid lines. The theorem of Kenyon and Schlenker Q asserts that for every such graph the covering lattice admits a rhombic 
embedding in the plane, that is one where all its edges have the same length. 

curves do not self-intersect and cross any other at most once) there is a rhombic embedding of C into the plane, that 
is, one in which all the edges have equal length. Each rhombus corresponds to an edge E of Q, and defines an opening 
angle as- As before, we may define parafermions on the edges of C and demand that they be discretely holomorphic 
when summed around the edges of each rhombus. This will be the case if we choose Xk{E) = Xkcicts), as given by 
(fT8)) . We conjecture that these values give the location of the FZ points on this inhomogeneous model. The rhombic 
embedding specifies how this model should be embedded in the plane in order that its scaling limit be given by the 
FZ conformal field theory. 

Finally we discuss how the conditions that allow the Yang-Baxter (star-triangle) relations in these models are 
compatible with discrete holomorphicity. Consider part of a lattice C where three lines meet at a point. This can be 
resolved in two different ways (see Fig.[5]). The Yang-Baxter relations guarantee that the Boltzmann weights, keeping 
all the other lines fixed, are independent of how this is done. If we 2-color the faces of C, it can be seen that the 
second resolution adds one more black face, that is one more vertex of Q. The invariance under this is the star-triangle 
relation. 

The rhombic embeddings in the two cases are shown in Fig. [51 It can be seen that in both cases the three rhombuses 
fit together to form a hexagon with opposite sides parallel, and that they simply correspond to two different tilings of 
the hexagon by the same three rhombi. The condition that the contour sum of the parafermions around each rhombus 
should vanish implies that the couplings should satisfy Xk — Xkci^^ — a), which are the critical star-triangle relations. 

IV. CONCLUSION AND DISCUSSION 

In this paper we considered the nearest neighbour Zjv model on some fairly general lattices. We identified the 
dual disorder operators, and used these to define parafermionic operators which reside on the edges of the covering 
lattice. For the isotropic square, triangular and honeycomb lattices we showed that, at the Fateev-Zamolodchikov 
critical points, these parafermions obey the discrete version of the Cauchy-Riemann relations, that is their contour 
sum around each elementary face of the covering lattice vanishes. We then extended this idea to regular lattices 
with anisotropic couplings, and showed that, if they are correctly embedded in the plane, it is possible once again to 
identify discretely holomorphic parafermions at particular values of the couplings. These we conjecture to correspond 
to the FZ points on these anisotropic lattice, which agrees with previously known cases. 

A crucial feature of these embeddings is that the faces of the covering lattice should be rhombi. This enabled us to 
extend our results to more general inhomogeneous lattices of the Baxter type, thanks to a theorem Q which guarantees 



8 




FIG. 5: Two different tilings of a hexagon by the same set of three rhombi. The right-hand case has, in the graph Q, an 
additional vertex associated to an order operator s as compared to that on the right. Discrete holomorphicity for each rhombus 
fixes the couplings on the dashed lines to be related by the star-triangle transformation. The two pictures are also related in 
the original graph £. by moving one of the curves past the vertex formed by the other two - the Yang-Baxter relation. 

the existence of a rhombic embedding in such cases. In this picture, the relation between discrete holomorphicity and 
the Yang-Baxter, or star-triangle, relations becomes particularly clear. 

We expect that the discretely holomorphic lattice parafermions become the holomorphic parafermions of the con- 
formal field theory, but this requires further assumptions. As pointed out in Refs. @, Q, if we regard the discrete 
Cauchy-Riemann equations as a linear system, there are in general twice as many unknowns as equations (only in the 
Ising case, when the phases of the parafermions are not free, are there the correct number). Morera's theorem, which 
allows one to deduce that a function / defined on is analytic if its contour integral around every closed contour 
vanishes, applies only if / is also assumed to be continuous. If we could prove that all correlators of our discrete 
parafermions become continuous functions in the scaling limit, this would be sufficient to show complex analyticity, 
at least with sufficiently smooth boundaries. However this step is highly non-trivial, as evidenced by the ambiguity 
in the interpretation of the linear relations discussed at the end of Sec. |TT1 There we showed that they lead to the 
identification of discrete parafermions iprf and also tprf — e^^^'^^'ij^rf. If they corresponded to conformal fields in the 
scaling limit, they would have conformal spins Pm and p,„ — 2. However, if the correlators of V' become continuous 
functions in the scaling limit, this cannot be true of those of if), and vice versa. 

If this problem can be overcome, we believe that our identification of suitable discretely holomorphic quantities 
should be the ffist step in showing that suitable defined curves in the Za? have SLE as their scaling limit, as has been 
recently conjectured [l3i[IH- In order to do this, following the ideas of Smirnov Q, however, it is necessary to identify 
these quantities with observables of these curves which are martingales of some discrete exploration process. Since 
the domain walls (or, equivalently, high-temperature graphs) of the Z^r model do not correspond directly to simple 
lattice curves (except for A'' = 2), there are a number of difficulties yet to be overcome in this program. 
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